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We study quantum critical behavior in three dimensional lattice Gross-Neveu models containing two massless 
Dirac fermions. We focus on two models with SU(2) flavor symmetry and either a Z2 or a U(l) chiral sym- 
metry. Both models could not be studied earlier due to sign problems. We use the fermion bag approach which 
is free of sign problems and compute critical exponents at the phase transitions. We estimate v = 0.83(1), 
r\ = 0.62(1), jty = 0.38(1) in the Z 2 and v = 0.849(8), r\ = 0.633(8), r^ = 0.373(3) in the 17(1) model. 
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The presence of massless Dirac fermions at low energies in 
graphene, has created much excitement over the past decade 
111 |2|]. By increasing the interaction strength between the 
electrons experimentally, an energy gap can be opened and 
fermions can become massive yl |4j]. Such quantum phase 
transitions, between a massless (semi-metal) and a massive 
(insulator) phase, are well known in particle physics. The 
possibility of studying them in a laboratory has ignited in- 
terest in the subject recently [5H7|]. Renormalization group 
arguments suggest that non-relativistic effects and long range 
interactions could be irrelevant [8], and the transition could 
belong to the universality class of similar phase transitions in 
three dimensional relativistic four- fermion field theories with 
two massless Dirac fermions [9, 10]. While Monte Carlo cal- 
culations of the critical exponents in models of graphene have 
emerged recently, the results are neither consistent with each 
other n 1 lMl 3fl nor do they match theoretical predictions [9]. 
Our current understanding of related four- fermion field the- 
ories is also quite limited. Compared to the precision with 
which three dimensional Ising and XY models have been stud- 
ied 1 14, 15], critical exponents in models with similar symme- 
try breaking patterns but in the presence of two massless Dirac 
fermions at the critical point remain largely unknown. As we 
explain below, some existing results are even puzzling. In this 
work we report new results in models that could not be studied 
earlier due to sign problems. Our results clarify some puzzles 
and help understand the subject better. 

Relativistic four-fermion models have a long history and 
are usually studied in the presence of either scalar interac- 
tions (Gross-Neveu models) or vector interactions (Third ng 
models) il6l421[l . Their lattice formulations using staggered 
fermions are popular, but due to fermion doubling one flavor 
of staggered fermions in three dimensions produces two fla- 
vors (Nf = 2) of Dirac fermions 11221 12311 . Symmetries of 
the microscopic models play an important role in determining 
the universality class of phase transitions. Gross-Neveu mod- 
els with a variety of symmetries have been studied using large 
Nf expansions 12412511 . e-expansions |26], renormalization- 
group (RG) flow methods ll27l - l29ll . and lattice Monte Carlo 
calculations IJ30H32I1 . Although much has been understood, 



there remain puzzles. For example, the critical exponents in 
the continuum Gross-Neveu model with a U(4) X Z^ sym- 
metry computed with the RG-fiow method 1 28], match those 
calculated with lattice Monte Carlo methods in a model with 
an SU(2) x Z2 symmetry IDOL Both models contain two 
flavors of Dirac fermions and calculations give v « 1.0 and 
i] « 0.75. Why do models with two different symmetries lead 
to the same critical behavior? Are symmetries dynamically 
enhanced in the lattice model at the critical point? On the 
other hand, is it possible that the results of Ref. |30J are incor- 
rect since sign problems were ignored 11331 13411 ? Here we show 
that another lattice model with the same symmetries give dif- 
ferent critical exponents. Another puzzle concerns a compar- 
ison between calculations of critical exponents in the contin- 
uum Thirring model with U (4) symmetry computed recently 
using the RG-flow method [35], and those in a lattice Thirring 
model with an SU(2) x 17(1) symmetry obtained with Monte 
Carlo calculations that do not suffer from sign problems [36- 
3811 . While both models contain two flavors of Dirac fermions, 
in the continuum one finds v w 2.4 and 77 « 1.4 while in 
the lattice one finds v w 0.85 and r\ ra 0.65. Although this 
disagreement can be attributed to the difference in the sym- 
metries, it does raise the question if lattice calculations have 
uncovered a new universality class. Here we show that lat- 



tice Gross-Neveu models defined in 12211 and lattice Thirring 
models defined in 113611 have the same symmetries and critical 
exponents. 

Lattice Gross-Neveu models with one flavor of staggered 
fermions cannot be solved reliably in the traditional approach 
due to sign problems J331I34I1 . The fermion bag approach is an 
alternative method which is free of sign problems and allows 
one to perform computations in these models reliably for the 
first time [38-40]. We use this new method to compute critical 
exponents in two types of lattice Gross-Neveu models invari- 
ant under either a Zi or a U(l) chiral symmetry. The models 
also have an additional SU(2) flavor symmetry which was 
appreciated only recently. Since they naturally describe two 
flavors of Dirac fermions in the critical region, these models 
have many properties similar to graphene including symme- 
tries. They were formulated originally with auxiliary fields 
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FIG. 1. A pictorial representation of the bond couplings Ul (left), 
Uw (center) and Um (right) discussed in the text. Each bond refers to 
the four-fermion interaction term of the form XxXx X y Xv 



that live at the center of cubes and couple to fermions on the 
corners [22]. After integrating over the auxiliary fields we ob- 
tain four-fermion models that couple fermion fields within a 
hypercube. Their action can be written as 



s = zZ^ x "> D ^v^y) 

%,y 



(xy) 



U{ X y)XxXx XyXy (1) 



where x(x),x(x) denote two Grassmann valued fermion 
fields at the lattice site x and D is the free massless staggered 
fermion matrix defined by 
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in which a labels the three directions 
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and r\x, a = 
= (1,1,0) are 



the staggered fermion phase factors 114 111 . The four-fermion 
interaction term involves the sum over three types bonds de- 
noted by (xy) (see Fig. [TJ: (1) link bonds L (where x, y are 
nearest neighbor sites), (2) face bonds F (where x, y are sites 
diagonally across faces of squares), (3) body bonds B (where 
x, y are sites diagonally across the bodies of cubes). 

In a general lattice four-fermion model the three couplings 
[/l, Uw an d Us will be arbitrary. However, in our study they 
are constrained since the action (Q3 is obtained by integrating 
over auxiliary fields from a model that contains a single cou- 
pling. In the Gross-Neveu model with Z2 chiral symmetry, 
we find £/l = 2L% = 4J7b = U, while with U(l) chiral sym- 
metry we find L?l = 4Um = U, Ur = [33]. In other words, 
face diagonal bonds break the 1/(1) symmetry to Z^. In ad- 
dition to chiral symmetries, models with action (fl3 have an 
SU(2) flavor symmetry which is hidden in the auxiliary field 
approach and was not appreciated earlier [42]. Indeed, when 
Up = it is easy to verify that the action (HJ is invariant under 
the following 5(7(2) x 17(1) symmetry, 



i;)-> e " y (£)' (*■ Xo)^(x Xo)Vh- i °, 

(3) 
where the subscripts e and o refer to even and odd sites and V 
is an SU(2) matrix. When (Jf^O the symmetry is restricted 
to 9 = ir/2 and the action is invariant only under an 517(2) x 
Z2 symmetry. 

Since four-fermion couplings are perturbatively irrelevant 
in three dimensions, models with action (Q]i have a massless 




FIG. 2. Plot of the chiral susceptibility at U — 00 for the Z2 (top) 
and U(l) (bottom) models. The solid curve in the top graph is a 
fit to the constant for L > 16, while in the bottom graph it is a fit 
to the finite size scaling form (J5) for L > 10 obtained from chiral 
perturbation theory. 



fermion phase at small couplings U. As the coupling in- 
creases, a second order phase transition to a massive fermion 
phase accompanied by spontaneous breaking of chiral sym- 
metries occurs at a critical coupling U c . Our goal is to study 
the critical exponents at this transition. However, before fo- 
cusing on the transition region, it is useful to understand qual- 
itatively the physics of the massive phase at large U. There is 
an important difference between spontaneous breaking of Z2 
and U(l) symmetries; the former does not produce massless 
Goldstone bosons while the latter does. It is important to dis- 
tinguish this feature in our results. For this purpose we have 
computed the chiral condensate susceptibility, 
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as a function of the lattice size L at U — 00. At infinite 
coupling our models can be mapped into a statistical model 
of closed packed dimers and can be updated efficiently us- 
ing worm algorithms [43]. Results obtained are shown in 
Fig. |2] As expected, finite size effects are enhanced in the 
[/(l) invariant model due to the presence of massless Gold- 
stone bosons. Results for L > 10 fit well to the leading order 
chiral perturbation theory form [44] 
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with S 2 = 0.844(1), p s = 0.381(3) and x 2 /d.o.f = 0.4. 
In contrast, the Z2 model shows very small finite size effects 
which indicates the absence of massless modes, and the data 
for L > 16 fits the constant 0.971(1) with a x 2 /d-o.f = 1.7. 
In order to uncover the properties of the quantum critical 
point we focus on the chiral susceptibility © and the fermion 
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0.0893(1) 0.83(1) 0.62(1) 0.38(1) 2.54(7) 9.33(5) 27.3(3) 55.3(1) 48.67(3) 34.4(1) -18.2(7) -51.2(6) 7.4(4) 259.2(10) 1.8 
0.1560(4) 0.82(2) 0.62(2) 0.37(1) 0.13(1) 0.09(1) 0.02(1) 0.004(1) 0.02(1) 34.0(1) -4.5(3) -1.4(3) -1.8(8) -0.5(2) 0.88 

TABLE I. Results of the combined fit of data in the critical region to Eqs. (0 in the Zi invariant model (top row) and (7(1) invariant model 
(bottom row). 




DC 



35 



30 



25 



- 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1— 




0.086 0.088 



0.09 
U 



0.092 0.094 



0.24 



0.18 



™_i 0.12 



0.09 



0.06 



- 1 — i — i — I — i — i — i — i — I — i — i — i — i — I — i — i — i — r - 




_i i i I i i i i I i i i i I i i i i_ 



0.145 0.15 



0.155 
U 



0.16 0.165 



+ 

CM 

of 




0.145 



0.165 



FIG. 3. Plots of xi ' L 2 v and RfL 2+v ^ as a function of U for L from 12 to 36. The solid lines show the combined fit which gives 
U c = 0.0893(1), v = 0.83(1), rj = 0.62(1) and ^ = 0.38(1) in the Z 2 case (top row) and U c = 0.1558(4), v = 0.82(2), 77 = 0.63(2), 
r]^ = 0.37(1) in the (7(1) case (bottom row). 



correlation function ratio 
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the following universal finite size scaling relations: 

4 
x/L 2 - v = y £fk[(U-U c )L$ 

fc=0 

4 
R f L 2 +^=Y,Pk[(U-U c )L± 
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where x is the origin or any translation of it by a multiple of where we have kept the first five terms in the Taylor series of 

two lattice spacings in each direction, and a is a unit vector the corresponding analytic functions. In order to compute the 

along each of the three directions. Since fermions are exactly critical exponents rj, v and r)^ we perform a single combined 

massless, in the vicinity of U c we expect \ and Rf to satisfy fit of the data in the critical region to Eqs. d7j with fourteen 



parameters. For the Z^ invariant model the combined fit of 
the data using lattice sizes ranging from 12 3 to 36 3 gives v — 
0.83(1), r) = 0.62(1), jfy = 0.38(1) and U c = 0.0893(1) 
with a x 2 /d.o.f. = 1.8. For the U(l) Gross-Neveu model, a 
similar combined fit in the same range of lattice sizes gives 
v = 0.82(2), r] = 0.62(2), r, 4 , = 0.37(1), U c = 0.1560(4) 
with a x 2 /d.o.f. — 0.88. Plots of our data along with the 
fits are shown in Fig. [3] The complete list of the fourteen 
fit parameters are listed in Tab. H From the results above, it 
seems like the critical exponents do not change much when 
chiral symmetries change from Z2 to U(l); the differences 
are small and lie within error bars. 

The critical exponents in the SU(2) x C/(l) symmetric lat- 
tice Gross-Neveu model obtained here, are also consistent 
with the exponents in the lattice Thirring model, which also 
has an action of the form (Q} except that C/l = U,Ug = 
Um = 113 811 . This is reassuring since the two models have 
the same lattice symmetries. Thus, calling one as the lattice 
Gross-Neveu model and the other as the lattice Thirring model 
is just a matter of taste. We can also study other SU(2) xU(l) 
invariant models by choosing a different set of couplings. Re- 
cently, the model with C/l = Um = U, C% = was also stud- 
ied and the critical exponents were again found to be similar 
04511 . Thus, it is tempting to combine all data from the three 
different studies and perform a single combined fit to extract 
a more accurate set of critical exponents. Using such a fit 
we estimate the critical exponents in the SU(2) x 17(1) sym- 
metric lattice models to be v = 0.849(8), r\ = 0.633(8) and 
th, = 0.373(3). 

Interestingly, the model studied in Ref. |30| is also an 
SU(2) x Z2 symmetric Gross-Neveu model. It is slightly dif- 
ferent from the model studied here since the auxiliary fields 
in the defining model live on sites instead of centers of hy- 
percubes. Integration over the auxiliary fields, which couple 
fermions on the six neighboring sites, produces four-fermion 
couplings of the form given in the action (fTJ with C/l = Us — 
0, Uw = U. However, in addition there is a non-zero next-to- 
nearest-neighbor four-fermion coupling along each direction, 
which is not present in our work. Since no lattice symmetries 
change, it seems very unlikely that these differences change 
the universality class of the phase transition. Hence, we be- 
lieve the critical exponents of the model studied in Ref. [30| 
should have been identical to our studies here. Unfortunately, 
this is not the case and we think that ignoring the sign problem 
in the auxiliary field approach could have distorted the results. 
It would be useful to repeat the calculation with the fermion 
bag approach. 

In this work we have been able to accurately compute the 
critical exponents at phase transitions in a class of 517(2) x Z^ 
and SU(2) x C/(l) symmetric four-fermion models involving 
two massless Dirac fermions in three dimensions. The critical 
exponents of the two models match within errors and a more 
accurate calculation is necessary to distinguish between them. 
Since the symmetries are different, we do not see any reason 
for the two exponents to be the same, however we are unable 
to rule out this possibility at the moment. As far as we can 



tell these critical exponents have not been verified in contin- 
uum field theory by the recently developed RG-flow method. 
However, we note that the e-expansion to second order in a 
Gross-Neveu model does agree with our results for the expo- 
nents v and r], but not for rj^ 12611 . Finally, given many sim- 
ilarity between graphene and staggered fermions, it would be 
interesting if the critical behavior in graphene falls in one of 
the universality classes studied here. 
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